Orbital ordering in Lao.sSri 5Mn04 studied by model Hartree-Fock calculation 
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We have investigated orbital ordering in the half-doped manganite Lao.5Sri.5Mn04, which displays 
spin, charge and orbital ordering, by means of unrestricted Hartree-Fock calculations on the multi- 
band p-d model. From recent experiment, it has become clear that Lao.5Sri.5Mn04 exhibits a cross- 
type {z"^ — x'^/y^ — z^) orbital ordering rather than the widely believed rod-type (3x^ — r"^ /^y'^ — r^) 
orbital ordering. The calculation reveals that cross- type (z'^ — x^ jy^ — z^) orbital ordering results 
from an effect of in-plane distortion as well as from the relatively long out-of-plane Mn-0 distance. 
For the "Mn^+" site, it is shown that the elongation along the c-axis of the MnOe octahedra leads 
to an anisotropic charge distribution rather than the isotropic one. 

PACS numbers: 71.30.+h, 75.30.-m, 71.27.+a, 71.20.Ps 



Perovskite-type transition-metal oxides have attracted 
much attention due to their wide variety of magnetic, 
electrical and structural properties. In particular, man- 
ganites are extensively studied because of their rich phys- 
ical properties such as colossal magnetoresistance (CMR) 
and spin, charge and orbital ordering 1]. Historically, 
their electronic states have been understood in terms 
of a double-exchange (DE) mechanism with localized 
t2g electrons and itinerant eg electrons |3|. Also, it 
has been recognized that the eg orbitals are strongly 
hybridized with the oxygen p orbitals, and participate 
in the cooperative Jahn- Teller (JT) distortion of the 
MnOe octahedra 0, [H . The degeneracy of the eg level 
is lifted by the JT distortion and the resulting orbital 
polarization at the Mn^+ sites can be accompanied by 
certain types of antiferromagnetic (AFM) spin order- 
ing. Half-doped manganites such as Lao.sCao.sMnOs and 
Lao.5Sri.5Mn04 exhibit a so-called CE-type AFM order- 
ing, as shown in Fig. 1. The magnetic moments of Mn 
on a zigzag chain are coupled ferromagnetically, whereas 
the neighboring chains are coupled antiferromagnetically. 
The charge and orbital ordering in the single-layered 
perovskite Lao.5Sri.5Mn04 and three-dimensional per- 
ovskite Ndo.sSro.sMnOs were investigated using resonant 
x-ray diffraction and an alternating Mn^+/Mn^+ pattern 
were observed directly '3,'^. Mizokawa and Fujimori 
showed that JT distortion consistent with orbital order- 
ing plays an important role in stabilizing the CE-type 
AFM state by means of unrestricted Hartree-Fock calcu- 
lations. The orbital ordering of the eg electrons had been 
thought of as a rod- type (3x^ — — r^) orbital order- 

ing, as shown in Fig. 1(a). More recently, however, from 
the studies of liner dichroism in the Mn 2p-edge x-ray ab- 
sorption 91 and resonant soft x-ray scattering at the Mn 
2p edge lOj of the Lao.5Sri.5Mn04, where Tco ^ 217 K 
and Tat 110 K 11], it was demonstrated that eg elec- 
trons exhibit predominantly cross- type (z^ — /y'^ — z^) 
orbital ordering as shown in Fig. 1(b). 

In order to understand the underlying mechanism of 
the cross-type orbital ordering, we have studied orbital 
ordering in Lao.5Sri.5Mn04 in the presence of lattice dis- 




FIG. 1: Models of orbital ordering in the Mn02 plane, (a) 
Rod-type (3x^ —r^/?>y^ —r^) orbital ordering; (b) Cross-type 



^ /y^ — z'^) orbital ordering ^9^. Also shown is the CE- 



type AFM ordering by arrows. 



tortion using unrestricted Hartree-Fock calculation on 
the Mn02 square lattice. We have employed the multi- 
band p-d model, where the tenfold degeneracy of the Mn 
3(i spin-orbitals and the sixfold degeneracy of the O 2p 
spin-orbitals are taken into account |l3, |l^ . 

The model Hamiltonian is given by 



H — Hp + Hd + H. 
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FIG. 2: (Color online) (a) Two parameters representing 
of lattice distortion at the Mn"^^, that is, the out-of-plane 
elongation of the MnOe octahedra represented hy R = 
c?out-of-piane/<iin-piane and the in-plane distortion of the Mn04 
plaquette represented by dL / din-p\ane- (b) Schematic picture 
of the lattice distortion in the Mn02 plane. 



_ .0 ^. 

i,a,'m>'m' ,a 



pd 



pd A 
klotra kotraa 



Pklo 



H.c. 



(3) 



(4) 



Here, d]^^^ is a creation operator for the m-th 3d elec- 
tron with spin a at site a of the ith unit cell and d\^^^ 
and p\i^ are creation operators for m-th 3d and /-th 2p 
Bloch electrons with spin a with wave vector /c, respec- 
tively. V^H, and 14^^^ are O 2p-0 2p and O 2p-Mn 3d 
transfer integrals. The transfer integrals are given in 
terms of Slater-Koster parameters (ppcr)^ {pp^)i {pda) 
and {pdn). In the present calculation, they are fixed at 
(ppa) = 0.60 eV, (ppn) = -0.15 eV, {pda) = -1.8 eV 
and (pdir) = 0.81 eV for the averaged Mn-0 distance 
of the Mn3+ site |i[ll[i3. Also, the atomic distance 
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assumed to obey Harrison's rule d~^-^ or d~'^ (ref. 'l^). 
The p-d overlap integrals are fixed at = 0.108 eV 
and aS'tt = —0.0486 eV. The intra-atomic Coulomb inter- 
action between the Mn 3d electrons is expressed using 
Kanamori parameters, u' and j ^3 5 which are fixed at 
u^v! ^ 2j, u = 7.19 eV, u' = 5.67 eV and j = 0.76 eV. 
The charge-transfer energy A is defined by — + 4/7 
for Mn^+, where and Cp are the energies of the bare 
3d and 2p orbitals and U{= u — 20/9j) is the multiplet- 
averaged d-d Coulomb interaction and is fixed at 4 eV. 
The unit cell for the CE-type AFM state thus contains 
eight Mn sites. The total electron number in the unit cell 
is 220 corresponding to the 8 Mn^-^+ and 32 ions. 

We have applied the Hartree-Fock approximation to 
the model Hamiltonian to simplify the Coulomb interac- 
tion terms to one-electron operators. For example, the 
second term in eq. (3) is expressed as 

^ ^ ^lam^^'ioim^dl^^^diami 
Q!,m 

i,a,m 

- «E *!^Ja™T^i""'T)(4amiC«i«mi>- (5) 

i,Q!,m 

The wave function is given by a single Slater determinant. 
We have iterated the self-consistency cycle until succes- 
sive differences of the order parameters (dl^^^diama) and 
{dlama^ic^m'a') Converged to less than 1x10"^. 

We consider two parameters for lattice distortion in 
Lao.5Sri.5Mn04, that is, the out-of-plane elongation ratio 
of the MnOe octahedra and the in-plane distortion of the 
Mn04 plaquette around the Mn^+ site as shown in Fig. 
2(a). The former is defined hy R = 6^out-of-piane/c^in-piane, 
where (iout-of-piane is the distance between Mn and out- 
of-plane oxygen and (iin-piane is that between Mn and in- 
plane oxygen. The latter is defined by (iL /c^in-piane, where 
dL {Mn^^) and ds {Mn^~^) are the longer and shorter 
Mn-0 distances respectively and satisfy the relationship 
2c^in-piane = c^L(Mn^+) + d5(Mn^+). The distortion in 
the Mn02 plane is shown in Fig. 2(b). 

In Fig. 3, we have plotted the distribution of electrons 
on the different Cg orbitals at the "Mn^+" site as a func- 
tion of the lattice distortion, R = C^out-of-plane/^^in-plane 

and (iL/<^in-piane- We have selected — y'^ and 3z'^ — 
as the basis orbitals for the Cg states. Although the lattice 
distortion causes a slight tilting of the MnOe octahedra, 
in the present model we have neglected the tilting and 
included only the bond length changes for simplicity. In 
the limit of wealc p-d hybridization, the electron num- 
bers of the Cg orbitals nx2_y2 = 0.75 and risz^-r^ = 0.25, 
that is, the relative population nx2_y2 /risz^-r'^ = 3 for 



the rod-type {3x 
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and risz^-r^ = 0.75, that is, nx2_y2 /n^^ 
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FIG. 5: (Color online) Difference in the number of electrons 
FIG. 3: (Color online) Orbital distribution of Cg electrons at between the "Mn^^" and "Mn'^^" sites 
the "Mn^+" site. 
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FIG. 4: (Color online) Orbital distribution of Cg electrons at 
the "Mn^+" site. 



for the cross- type {z^ — /y^ — z^) Therefore, 
one can use nx2_y2 /n^^2_^2 as a measure of the type 
of orbital ordering. When the in-plane distortion is 
absent (<iL/<^in-piane = 1) and the ratio of the out-of- 
plane to in-plane Mn-0 distances is small {R ^ 0.9), 
the — orbital is dominantly populated. As long 
as the in-plane distortion at the "Mn^+" site is weak 
(o^L/c^in-piane — 1) and R < 1.0, the orbital distribution 
of Cg electrons remains close to the rod-type orbital or- 
dering. However, when the in-plane distortion becomes 
strong (<iL/<^in-piane > 1)? the distribution of Cg electrons 



approaches that of the cross-type orbital ordering. With 
the in-plane distortion and the elongation along the c- 
axis, Lao.5Sri.5Mn04 shows stronger tendency toward an 
orbital polarization of strong z character. The dominant 
occupation of the 3z^ — orbital causes the abrupt drop 
as shown in Fig. 5. In the previous HF calculation, which 
includes only the in-plane distortion, it has been reported 
that the "Mn^+" sites in the CE-type AFM states are 
accompanied by the rod- type orbital ordering . In the 
present case, by considering the out-of-plane elongation 
in addition to the in-plane distortion, we have shown 
that both types of distortions are important to realize 
the cross-type orbital ordering. We have also carried out 
calculations by changing the Coulomb interaction [T^ . 
As a result, the relative population of the eg orbital for 
"Mn^+" changes slightly. However, the orbital ordering 
has also been found to be mainly determined by the ef- 
fect of lattice distortion rather than that of the on-site 
Coulomb interaction. 

As for the three-dimensional Mn perovskites, the res- 
onant x-ray scattering study of Ndo.sSro.sMnOa has 
strongly suggested an orbital ordering of the rod-type 
(3x^ — r^/37/^ — r^) 0. Because the Mn-0 distance along 
the c-axis in the Ndo.sSro.sMnOa is much shorter than 
that within ah plane, the — y'^ orbital becomes occu- 
pied by more electrons, and the rod-type orbital order- 
ing is favored. This is contrasted with the present case 
of Lao.5Sri.5Mn04, where the in-plane and out-of-plane 
Mn-0 distances are similarly strong, and the cross-type 
orbital ordering is realized. 

We have also calculated the orbital distribution of eg 
electrons at the "Mn^+" site of Lao.5Sri.5Mn04. In Fig. 
4, we have plotted the orbital distribution of eg electrons 
at the "Mn^+" site as a function of lattice distortion. 
Here, it should be noted that there is finite number of 



4 



Cg electrons at the "Mn^+" site in spite of the formal 
configuration because of the p-d hybridization. For 
R < 1, the orbital distribution at the "Mn^+" site is 
isotropic {71^2 _y2 :^ 713^2 _^2) as expected for the con- 
figuration of the "Mn^+" ion. However, for > 1, the 
— orbital becomes partially occupied through the 
anisotropic p-d hybridization of the distortion MnOe oc- 
tahedra, and the orbital distribution at the "Mn^+" site 
becomes anisotropic. In fact, because the MnOe octahe- 
dra for the "Mn^+" site in Lao.5Sri.5Mn04 is elongated 
along the c-axis, the Mn-0 distances within the ah plane 
become short in comparison with that along the c-axis. 
As a result, the strong in-plane covalent bonding leads to 
a net charge transfer to the — y'^ orbital. As R further 
increases, 71^2 _y2 strongly increases, as indicated by the 
sudden increase of 71^2 _y2 / 712, z'^ _j.2 in Fig. 4. Because 
of the p-d covalence, the net charge difference between 
"Mn^+" and "Mn^+" is significantly reduced from 1 as 
shown in Fig. 5. For large i?, the difference in the num- 
ber of electrons between "Mn'^+" and "Mn^+" sites is 
reversed as shown by the sudden drop in Fig. 5. This 



means that the CE-type AFM state in Lao.5Sri.5Mn04 
becomes unstable for strongly elongating of the MnOe 
octadedra along the c-axis. 

In conclusion, we have studied the orbital ordering in 
Lao.5Sri.5Mn04 by means of unrestricted Hartree-Fock 
calculations. It is shown that cross- type {z^ —x^ jy^ — z^) 
orbital ordering originates from the in-plane distortion of 
the MnOe octahedra as well as the relatively long out-of- 
plane Mn-0 distance. Further, we have found that the 
orbital distribution at the "Mn^+" site is isotropic for 

< 1, but, becomes anisotropic for > 1. Other phys- 
ical implications of the different kinds of distortions and 
the orbital ordering in manganites are an open question 
for future studies. 
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